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Using the general framework of quantum field theory, we derive basic equations of quantum field 
kinetics. The main goal of this approach is to compute the observables associated with a quark- 
gluon plasma at different stages of its evolution. We start by rewriting the integral equations for the 
field correlators in different forms, depending on the relevant dynamical features at each different 
stage. Next, two versions of perturbation expansion are considered. The first is best suited for 
the calculation of electromagnetic emission from chaotic, but not equilibrated, strongly interacting 
matter. The second version allows one to derive evolution equations, which are generalizations of 
the familiar QCD evolution equations, and provide a basis for the calculation of the initial quark 
and gluon distributions after the first hard interaction of the heavy ions. 
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O ' 1. Introduction 

^ The most ambitious goal of the RHIC and LHC programs is to discover a new state of the matter - the quark-gluon 

plasma (QGP). Evidence for a QGP will require a self-consistent analysis of many signals from all stages of the collision 
■ - initial "hard" processes (r < 1/2 fm/c), the QGP itself, and from the "cool" hadronic gas (r ~ 10 — 30 fm/c). 
[ Thus, a continuous description for all t is of more than of academic interest. This is a difficult task, and currently, 
each stage is described using different approach. Here, we primarily wish to design a formalism that allows one to 
^-H ' describe all stages of the collision, including the transient ones, using the same technical tools. 

K*" . The existence of a QGP is inseparable from the process of creation of the matter it consists of. It is an essentially 
' quantum process . Therefore an exact definition of the initial state of the system, and of the observables in the 
, expected final state, is required. Two difficulties arise: (i) It is unclear how the stable nuclei of the initial state are 
' build up from quarks and gluons; and (ii) the expected final state is imperfectly understood. Thus, the theory should 
' be adaptable enough to deal with these uncertainties. It must bridge the logical gap between the language used to 
t describe stable nuclei and a QGP, as is done phenomenologically in using an intermediate free parton language. More 
On ' specifically, in describing a transition between the initial state of two stable, well shaped nuclei, and a system of free 
, quarks and gluons, we encounter one of the most painful problems of quantum field theory - initial and final states 
O i' of the system can not be described in terms of the same language. Before the initial collision, the system is confined 
' , and the vacuum is dominated by quark and gluon condensates, while afterwards all the condensates are destroyed 
and the quark-gluon dynamics is calculated with respect to the perturbative vacuum. This vacuum itself is a product 
(— I I of the collision. 

• • • We may begin the study with the assumption that a high multiplicity of quarks and gluons has already been created. 
. 5^ ' However, this "plasma" may remain out of thermal and chemical equilibrium for a long time, and may not reach it 
. at all In this case, the calculation of even relatively simple signatures such as dilepton and photon emission is a 
' nontrivial problem. An overview of the many publications on this subject reveals that the QGP has generally been 
5^ , treated as a totally equilibrated system, and that calculations are essentially based on the detail balance relationships. 

Therefore, we may conclude that any reasonable theory should rely neither on detail balance and thermal equilib- 
rium, nor even on the existence of the same ground state for all stages of the QGP evolution. The theory should 
explicitly follow the temporal sequence of the stages, and allow for smooth transitions between them. The prototype 
of such a theory was designed by Keldysh j^] for non-equilibrium condensed matter systems with the aim of deriving 
a quasi-classical kinetic equation. 

The version of the Keldysh technique presented below, and named "Quantum Field Kinetics" (QFK), takes an 
intermediate position between the theory of scattering and quantum field theory of the many-particle systems. It 
contains both as its limits. On the one hand, it allows one to calculate the inclusive cross-sections and rates of 
emission since the summation over all unobserved states is implicit. Also, it provides the proper balance between 
"real processes" and "radiative corrections" which results in the cancellation of infrared divergencies, and makes 
unnecessary to introduce artificial intermediate cut-offs. On the other hand, QFK allows one to calculate local 
observables for many-particle systems. As a by-product it recovers microcausality, precisely in the form which was 
used in axiomatic field theory. 

For any particular problem, QFK-based calculations always begin with the basic definition of observables in terms of 
their Heisenberg operators, the Lagrangian of the theory, and the density matrix of the initial state. Several examples, 
motivated by the expected scenario in heavy ion collisions and possible probes of the QGP, are described in Section 
2. One may express the observables via field correlators with differently ordered field operators. The correlators are 
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naturally arrayed in 2 x 2 matrices, and we follow the original Keldysh idea of contour ordering. The correlators 
obey a system of matrix equations of the Schwinger-Dyson type, which is easily derived via functional methods. The 
equations are reviewed in Sec. 3. 3. 

Starting with the initial form of the matrix equations, it is straightforward to generate a formal perturbation 
series for the probabilities of inclusive processes (instead of amplitubes for exclusive ones) in powers of the coupling 
constant. This series reproduces the Feynman expansion diagram for diagram, except that now each vertex acquires an 
additional dichotomic index. The way to sum this perturbation series is always influenced by the physically motivated 
renormalization conditions, which may be different from case to case. 

As a preliminary step for the future rearrangement of the perturbation series, we begin, also following Keldysh, with 
a rotation of the matrix basis. We introduce the retarded and advanced Green functions of QCD along with additional 
correlators which carry information about the space of states and derive the equations they obey. In Section 3.4 we 
show that the latter allow for a formal solution which comes to be the first approximation of the QCD evolution 
equations, if we renormalize the retarded propagators according to the requirement of light cone propagation. 

It is important to note that instead of using the traditional method of the operator product expansion, we can 
derive the evolution equations for the DIS structure functions using the language of the QFK. This language avoids the 
parton phenomenology which is needed when DIS structure functions are used for the computation of the quark and 
gluon distributions in first hard collision of the two nuclei. Thus, QFK provides a firmer footing for the description 
of heavy ion collisions. 

In Section 4 we discuss different versions of the perturbation expansion. We show that the ordinary perturbation 
series in powers of the coupling constant is adequate for calculations of systems with slowly varying macroscopic 
parameters. We describe the scheme using the example of dilepton emission from the QGP. An expansion which 
is suitable for the violent impact of relativistic composite systems requires another kind of expansion, one which 
preserves the leading light cone singularity of the propagators. 

2. Signals from the quark-gluon plasma 

In this section we begin the design of a technique which will allow us to calculate different observables associated 
with QGP. The latter should be defined unambiguously at both the theoretical and apparatus level. We shall try to 
trace their origin using only the basic principles of quantum mechanics for as long as is possible. If successful, we may 
then express the signals in terms of the parameters of the emitting system. All the systems we shall study are described 
by the standard QCD Lagrangian. In order to deal with electromagnetic probes we include the interaction of charged 
quarks, q{x), and leptons, ipix), with the photon field, yl'^(x), in the total renormalized interaction Lagrangian: 

q,i q.i 

+grfabcd^Bl{x)Bl{x)Bl{x) + {9llA)}abcfaghB^,(x)B':^{x)Bl(x)Bt{x), (2.1) 

where B^(x) is the gluon field, and gr is a renormalized strong coupling constant. Other notation is commonly used, 
and requires no comment. The standard counterterms which will be used for renormalization are as follows: 

>Cct(x) = ^^Bi^{g^,d' - d^d, + Xn^n,)B: - {Z, - l)grfabcd^B:{x)Bl{x)Bt{x) - 

^{Zi - l)^fabcfagkB^(.x)B'^{x)Bf^{x)Bl:ix) + [Z^ - l)Ux)irdM^) + {ZiF ~ l)gMx)t%r B';{x)qj{x) (2.2) 

Here, all {Z — l)'s are considered as the small parameters. In what follows we shall either use covariant gauges in the 
lowest order calculations when the ghosts do not contribute yet, or shall work in the axial gauge, B'^Uf^ = 0, where 
they are absent. 



The Lagrangian ( p.l| )- (2.2) gives rise to the ordinary S-matrix in the m-interaction picture, 

S = TeMi J d^xC„,tix)} (2.3) 

which is considered to be a limit of the evolution operator that governs the dynamics of Heisenberg observables. 

The initial state for any system can be described by the density matrix pqcd- It is formed due to strong interactions 
only, and we shall specify it later. When the inclusive cross-sections of photon and dilepton emission are chosen as 
the observables, this density operator should be augmented by the projector on the initial vacuum state of photons 
and leptons, 
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P = PQCD ® |0e7)(0eT,|, 



(2.4) 



and we assume that the QGP remains transparent for photons and leptons throughout its history. The amphtudes of 
the transition from one initial state, |m), to a final one containing photon or dilepton read as 

(X|c(k,A)S'|m) or {X\h{2)a{l)S\in) , (2.5) 

where c(k, A) and a( J) — a(kj, uj) and 6( J) are the photon, electron and positron annihilation operators. Summing 
the squared moduli of these amplitudes over a complete set of uncontrolled states and averaging over the initial 
ensemble, we find the inclusive spectra of photons and dileptons, 

^^Y. Spp.n^^ct(k, A)c(k, \)S, (2.6) 

A 



(ikidk2 



(2.7) 



It is easy now to commute the Fock operators with S and S'^ . Due to the QED vacuum projector in the density 
matrix (2.4), only the commutators survive and Eq.(|2.6|) takes the form; 



j3 dN, 



dkd^x 2{2n 



^<(-fc), 



where k'^ > 0, and the polarization tensor 



5S 



M^(x) SA-'iy)' ' 



(2.8) 



(2.9) 



is the Fourier transform of the product of two Heisenberg electromagnetic currents averaged with the density matrix 

PQCD- 



SS 



(2.10) 



The Heisenberg current (as any other local Heisenberg operator) can be written down in the two equivalent ways, 

r(a;) = S^T{r{x)S) ^ T^{r{x)S^)S, (2.11) 

where j^{x) = (1/2)^^ ^eq[qi{x)j^,qi{x)] , is the operator of the electromagnetic current in the m-interaction 
picture. Relations such as (2.11) are extremely important as they allow one to keep the initial order of the operators 
through all stages of calculations. This order is strictly prescribed by the definition of the observables and may not 
be changed safely, except under very spe cial circumstances. 
Now the dilepton rate of emission (2^) takes the form 



^ ^ dkidk2C?'*a; 



. 2 Lfj_,j(ki , k2) ^ 
4(2^)6 



A5'o(-fc), 



(2.12) 



where fc = fci + ^2, and L^^'^ — kf^k!^ + kl^k^ — g'"'(fcifc2 — m^) is the polarization sum of the lepton spinors. The 
electromagnetic correlator 

A^;(-fc) = -* / dHx-y){T\A>^{x)S^)T{A>^{y)S))e-'^''^'^-y\ 



(2.13) 

is a kind of photon Wightman function averaged over p^'^^ . The opera tor A{x) of the in- interaction picture and 
the Heisenberg operator A(a;) are connected via relations similar to (2.11). In the absence of radiative corrections to 
the photon propagation, dyna mical equations which will be derived in the next Section will allow us to rewrite the 
photon correlator in Eq. ( ^.12 ) as following. 



(2.14) 
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In order to derive Eqs. (2.S) and (2. IS) we have assumed that an exphcit separation of long-range and short-range 
scales is possible, and have introduced the emission rates per unit volume, instead of the inclusive cross-sections. 

The way we proceeded above demonstrates the very simple and general principle of formulation of the problem: 
Obscrvables like cross-sections or the rates of emission must be expressed in terms of specifically ordered products of 



the Heisenberg operators. If we change the photon and Icpton creation and annihilation operators in the Eq. (2.5) 



( p. 7]) for those of quarks or gluons, and replace the pi„ of the QGP by a density matrix for the two colliding nuclei, 
we obtain the starting point for a computation of the quark and gluon distributions after the first hard collision. 

We now require a formalism which will allow us to calculate these quantities. We are to keep in mind that all 
the operators are driven by Cmt, and that all the information about the initial state of the system is hidden in pi„. 
Nothing else is needed to solve the problem. 

3. The equations of relativistic quantum field kinetics 



3.1. Basic definitions 

In this study, the calculation of obscrvables such as the emission rate is based on work by Keldysh ^ . It incorporates 
a specific set of exact (dressed) field correlators. These correlators are products of Heisenberg operators, averaged 
with the density matrix of the initial state. For the quark field they read 

Gio(a:, y) = -i(q(a;)q(y)), Goi(a;, y)= i(q(2/)q(a;)), 

Goo(2:,y) --z(r(q(a:)q(y))),Gn(a:,2;) = -*(Tt(q(a:)q(2/))), (3.1) 

where T and are the symbols of the time and anti-time ordering. They may be rewritten in a unified form, 

GAB(a;,y) = -^(T,(q(xA)q(2/B))) , (3.2) 

in terms of a special ordering Tc along a contour C = Co + Ci, the doubled time axis, with T-ordering on Co and 
Tt-ordering on Ci. The operators labelled by '1' are T^-ordered, and stand before the T-ordered operators labelled 
by '0'. RecaUing that 

q(x) = S^T{q{x)S) ^ T\q{x)S)S, (3.3) 

we may introduce the formal operator Sc — S, and rewrite ( |3.2| ) using the operators of the m-interaction picture 

GAB{x,y) = -i{T,{q{xA)q{yB)Sc)), (3.4) 

where, by the definition, the internal variables of S lie on Co, and those of on Ci. 

The boson Greenians are built in the same manner, i.e., for gluon field, B{x), and photon field, A(x), we have 

T>ab{x, y) = ~i{T,{B{xA)B{yB)S,)), (3.5) 
Aab{.x, y) = -i{T,{A{xA)A{yB)Sc)), (3.6) 

where the vector and colour indices have been suppressed. 

We do not consider the path C to be extended to complex values of the time i, nor to be closed. Moreover, to some 
extent we take a step backwards by restoring many elements of the old-fashioned perturbation theory which deals 
with retarded and advanced propagators, and incorporates microcausality in the sense of Yang-Feldman equations 

. The C-contour technique is only a convenient trick to do this in an economic way. Similarly, the T-ordered Green 
function is the sole carrier of physical information only under very special circumstances. Indeed, by definition 

Goo{x,y)^~iTY[pS^T{q{x)q(y)S)] . (3.7) 

If the density matrix of the system, p, corresponds to the exact stationary ground state, then in presence of an 
interaction the state-vector can acquire only a phase factor, and we obtain Feynman's Green function: 

I ^ -i^^pT{q{x)q{y)S) 

GF[x,y) = (3.8) 

(o)o 

As we shall see in a while, the Green functions Goo(a^, y) and Gpix, y) even obey different integral equations. 
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3.2. The density matrix of the initial state 



The choice of the density matrix specifies the physical phenomenon under consideration. If the initial state of the 
system consists of a only few excitations of the vacuum, then we have the density matrix of a pure state. In this case 
we are dealing with the well known picture of scattering. Eventually, this situation is described by a certain set of 
vacuum expectation values, and the corresponding density matrix of the pure ground state is 

PQCD = \0qcd){0qcd\- (3.9) 

The various density matrices which we shall use later for the computation of quark and gluon production in the first 
hard AA-coUision, and for the rates of the photon and dilepton emission from the QGP, emerge from the following 
scenario for the heavy ion collision: In the initial stages of a collision at RHIC or LHC energies (t ~ O.lfm/c), hard 
impacts take place. The initial state that precedes these processes consists of two stable, well shaped nuclei moving 
along opposite directions of the light cone, and the density matrix for each of them is the same as in the deep inelastic 
e-p or /LtA scattering. At this stage, the most interesting observables are the quark and gluon distributions after 
the first interaction that destroyed the nucleons. The QGP domain begins somewhat later when the distributions 
are already chaotized, and may be described by one-particle distributions. The most general density matrix which 
simulates any given form of a one-particle distribution is of the following form: 

p = nn^^^''*'^'''^"^'''^''^"''^'^''^ (3.10) 

N p,j 

where N labels the space cells at the hypersurface of the initial data, and nj{N,p) = aj{N,p)aj{N,p) is an operator 
of the number of partons of type j and quantum number p in the A^-th cell. Thus, we completely neglect all correlation 
effects in the initial phase space. Introduction of the cells is necessary only for extended objects without long-scale 
quantum coherence. For a single hadron they are not needed. 

The density matrix of the Gibbs ensemble of noninteracting quarks and gluons against a hydrodynamic background 
is of the same kind as ( 3.1C| ), and obeys the additional condition of being an entropy extremum. Introducing the local 



4- velocity of continuous media u'^{x) we may write 

P CD - TT "^^Pt^"^^"^ + PnQn)/Tn] (3 11) 

^'^'-^^ ^J- Sp[exp[{~PNUN + fJ^NQN)/TN]] 

where and Qat are the total 4- momentum and (baryonic) charge of free quarks and gluons at temperature T/v, 
and /iAf is the chemical potential in a small 3-volume Vjv on the hypersurface of the initial data. 

The explicit form of the "bare" Greenians which are the basis of the perturbation theory is quite evident: The 
"vacuum" Green functions are of the standard form, 

G%\,is)^ -2mi/^ + m)ei±so)6is^-m^), G*"^^ (s) = ±^-^±^ , 

DZoii^) = - 2md^'^is)ei±so)S{s^) ,<5'r(^) = ■ (3-12) 

where the projector d^'^{s) depends upon the choice of gauge. Usually, the bare Greenians of the ensemble are of the 
form 

Gab{s) = d^l{s) + Gp{s), D^/^is) ^ Dfri^) + ^r(^)' (3-13) 

where the additional terms originating from the 

Gp{s) = 2m{/k + m)6{s'^ - m2)[6'(so)nW(s) + 0(-so)n(-)(s)], 

Df{s) = -27rzd^''(s)J(s2)[0(so)/W(s) + 0(-so)/(-)(s)], (3.14) 

manifestly contain the Fermi- or Bose-occupation numbers, n^^^ or Z^^-*, respectively. They are the diagonal elements 
of the density matrix. All theorems of the Wick-type, which are necessary for the calculations, can be proven easily 



for density matrices like (3.10) and (3.11) 



3.3. The Schwinger-Dyson equations for QFK 
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Except for the matrix form, the Schwinger-Dyson equations for the Heisenberg correlators remain the same as in 
any other technique. An elegant and universal way to derive them which does not rely on the initial diagram expansion 
can be found in Ref . . For the quark field these equations are of the form 

Gab = Gab + ^ Gar ° ^bs ° Gss, (3.15) 

RS 

Here the dot stands for convolution in coordinate space, and for the usual product in momentum space (providing the 
system can be treated as homogeneous in space and time). Indeed, the only tool used to derive these equations was the 
Wick theorem for the ordered products of the operators. The type of ordering is inessential Explicit expressions 
for the self-energies are obtained automatically in course of the derivation of the Schwinger-Dyson equations. The 
quark self-energy matrix reads 

^AB{x,y) = i(-l)^+^g.' J2 (-1)''+^ / didrit^rGARix,OrjiB,sii,y;v)'DfA,Xf.iv,x) . (3.16) 
The form the strong qqB-vertex that appears is 



r^'<^,.(^,y;^)-(-ir^'+^ ^ ^ , (3.17) 



which is the functional derivative of the inverse Greenian of the quark field with respect to the "external" gluon field 
The matrix of the photon correlators obeys similar equations, 

Aab ~ ^AB + J2^aroIIrs° AsB, (3.18) 

RS 

where the electromagnetic polarization operator is 

7'^'^(x,y) = -z(-l)^+V' E i-lf^' fd^dr^rGAR{x,OE-ks,B{^'V;y)GsA{v,x), (3.19) 

and where the electromagnetic vertex 

E^BsA-^ y; .) = i-i)n.s,R S[G-Hy)]n, ^^^^^^ 

edAx[zp) 

is dressed by the strong interaction. The latter in turn obeys the equation 

= 7^ + £;^oGGoK, (3.21) 

with a four-fermion vertex K. 

The gluon field correlators also obey similar equations 

T>AB^ Dab + Y.^aR°^RS°'Dsb, (3.22) 

RS 

where the gluon self-energy has two terms: 

1 



^'ATix,y) = -igl E (-l)^+^+^+^[ / d£,dr^TTrt'^GAR{x.Oy'^s,B{^^V:v)GsA{v,x) - 

R,S=0 •' 

- J d^driV:^:f\x,^v'W^'R''^{i,a^''^^^^^^^^^ . (3.23) 



The 3-gluon vertex (in coordinate representation) is defined as 

[] 



n^lRSpi-:y^^) = (-l)^+-+^^PlglMf^ . (3.24) 
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We have omitted the trivial term with the 4-gluon vertex in Eq. (3.23). The exphcit coordinate expression for the 
bare 3-gluon vertex is 



fcc;l//3 



gr6Bi{zp) 



(3.25) 



where the inverse Green function ^(q) denotes a hnearized differential operator of the wave equation for the gluon 
in an external gluon field. The only importance of 3-gluon vertex in this form is that it displays its local nature, 
Vabc{x, y, z) dAB^AC^ix — y)S{x — z). In momentum representation it may be written as 



^ABC;a6cbl,P2,P3) = -*gr<5AS<5Ac/'^'l5"^(pi " Ps)^ + /^(p2 " Pa)" + ff^"(p3 " Pi)''] 



(3.26) 



The four types of operator ordering which enter Eqs. (3.1) are not linearly independent, i.e., there exists a set of 
relations between the field correlators, and between the self-energies: 



Goo + Gil — Gio + Goi, Soo + ^ii — ^^lo ^ ^oi, 
^00 + Dii — Dio + Doi, IIoo -I- IIii — — IIio — Hoi 



(3.27) 



These indicate that only three elements of the 2x2 matrices G, S, etc. are independent. To remove the overdeter- 
mination let us introduce new functions 



Gret — Goo ^ Goi, Gadv — Gqq — Gio, Gi — Goo + Gn; 
Sret = Soo + Soi, ^adv — ^00 + ^lo, Si ~ Soo + Sn, 



(3.28) 



as well as their analogs for bosonic correlators. One of the possible ways to exclude the extraneous quantities is to 
use the following unitary transformation |^ , 



G — R ^GR, S — R ^Si?, R — — ^ 

%/2 



1 1 

-1 1 



In this new representation, the matrices of the field correlators and self-energies have a triangle form, 

G 






Gadv 


, Af = 


Si 


^ret 


Gret 


Gi 








(3.29) 



(3.30) 



Applying transformation (3.29) to the matrix Schwinger-Dyson equations (3.15) we may rewrite them in the following 
form: 



Gi 



Gret — Gret + Gret ° ^ret ° G^-et, 
Godu = Gadv + Gadv ° ^adv ° Gadu; 
Gi + Gret ° ^ret ° Gi + Gi O Yiadv ° Gadv + Gret O Si O Gadv 



(3.31) 
(3.32) 
(3.33) 



There are now two ways to proceed with a further study of Eq. (3.33). If, as typical of condensed matter systems, 
one has two well separated short- and long-range scales, then it is reasonable to look for a kind of differential kinetic 
equation 0. The time- and space-derivatives on its l.h.s. correspond to the slow variations of collective parameters, 
while the short-distance dynamics is absorbed into the collision term. The existence of two separate scales here 
is considered part of the external physical input. It must be either self-evident, or be proven by a separate study 
(confirmed by observations). In order to derive the quasi-classical kinetic equation, one acts on Eq. (3.33) separately 
from the left and from the right with the differential operator of the free wave equation, and takes the difference of 
the two resulting equations. Further, it is useful to transform Gi to Wigner variables, which effectively separate the 
short- and long-range scales. This procedure results in an integro-differential kinetic equation for the density Gi in 
phase space. It is differential with respect to the long scale, and integral with respect to the short scale. 

In relativistic Boltzmann-type kinetic theories one encounters specific difficulties: First, the phase-space Wigner 
distributions are no longer positive-definite, and thus do not carry any direct physical information. Second, for 
phenomena like heavy ion collisions, the two scale dynamics is not evident a priori. Thus, it is safer to accept the 
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integral equations (3.31)-( 3.33| ) as they are; however, in th is forrn they are badly suited for practical calculations. In 
a short while we shall show that the exact equations ( 3.31 )^( ^.33 ) can be solved (at least formally). 

3.4- The formal solution of the integral equations 

The solution we shall look for now is equivalent to the rearrangement of the perturbation series for observables. It 
is useful when specific features of the non-equilibrium system must be taken into account. First, let us introduce two 
new correlators: 

Go — Gret — Gadv = Gif) — Gqi, (3.34) 

which coincide with the anti-commutator of the quark fields and thus disappear outside the light cone, and 

Eg = S,-et — '^adv — " ^iQ — SqI- (3.35) 



which is the commutator of two fermion sources an d has the s ame c ausal properties as (3.34). The integral equation 
for Go may be derived by taking difference of Eqs. (3.31) and (3.33): 



Gq + Gret ° ^ret ° Gg + Gq O T,adv ° Gadv + Gret O Eq O Gq, 



dv • 



(3.36) 



The sum and the difference of Eqs. ( 3.33 ) and ( 3.36 ) give corresponding equations for the off-diagonal correlators Gio 
and Gqi: 



Gni — Gni + Gret ° ^ret ^ 



Goi O ^adv ° Gadv — Gret ° 5^01 O Ga 



dv • 



Since Eq. ( 3.31 ) for the retarded propagator may be identically rewritten in the same form, 

Gret — Gret + Gret ° '^adv ° Gadv + Gret ° '^ret ° Gret — Gret ° '^adv ° Gadv, 

we may use Eq. ( 3.37| ) and derive the corresponding equations for the T- and T^-ordered propagators: 

Goo = Goo + Gret ° ^ret ° Goo + Goo o ^adv ° Gadv + Gret o En o Gadii- 



(3.37) 



(3.38) 



(3.39) 



This chain of routine transformations reduces equations for all elements of the matrix correlator Gab to a unified 
form. On the one hand, this representation shows that linear relations between the correlators (or, explicitly, different 
types of orderings) hold even for the equations that the correlators obey. On the other hand, this representation of the 
equations singles out the role of retarded and advanced propagators over all other corre lator s. In order to understand 
why their role is special, let us transform them further, and begin by rewriting of Eq. ( 3.33 ) for the density of states 
Gi identically as: 



(1 — Gret ° ^ret) O Gi — Gi O (1 + T^adv ° Gadv) + Gret O Ei O Gadv ■ 

Since Gret o ^ret ° Gret = Gret — Gret , it is casy to show that 

(1 + Gret ° 5]ret)(l — Gret ° '^ret) — 1 

Further, we have the following two relations 



(1 + ^adv ° Gadv) ^^(o) °Gadv , (1 + G^ef O E^gf) — Gret° Gj^ , 



(0) 



where 



^(0) (^) />x-m, Gf^] (x) = -i 



(0) 



(3.40) 
(3.41) 

(3.42) 

(3.43) 



are the left and right differential operators of the Dirac equation, respectively. Now, multiplying Eq. ( |3.40| ) by 
(1 + Gret o Sret) from the left, we find the final form of the equation we are looking for: 



Gi — Gret° Gj.g^ oGiO G^p"J °Gadv + Gret O Ei O Gadv 



(3.44) 
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Repeating these transformations for the other equations, we obtain the corresponding forms that are most convenient 
for subsequent analysis: 



Gio - Gref ° ^{0) 



°Gadv ^ Gret ° Sio O Gq 



dv : 



(3.45) 



Goo ~ Gret^ ^(0) ^(0) 



Sii O Ga 



dv • 



(3.46) 



While equations ( ^.37] ) and (3.39) are standard integral equations which have the unknown function on both sides, 
after transformation the unknown function appears only on the l.h.s., and we may consider Eqs. ( |3.44 )-(3.46) as the 
formal representation of the required solution. 

At this point, the first and the most naive idea is to ignore the arrows indicating the direction the differential 
opera tors act in, and to rewrite Eqs. ( |3.45| )-( 3.46 ) in momentum representation. Then the first term in each of 
Eqs. ( 3.45| ) will contain the expression {p^ — mF)5{p'^ — m?), which equals zero. It reflects simple fact that the off- 
diagonal correlators Gio and Gqi are solutions of the homogeneous Dirac eq uatio n. Ho wever , this approach does not 
appear to be sufficiently consistent: we would lose the identity between Eqs. (3.45) and (3.37), and make it impossible 
to generate the standard perturbative expansion in powers of the coupling constant. Eqs. ( 3.46 ) will be corrupted 
also. 



A more careful examination of Eqs. ( 3.45 )-( 3.46 ) shows that all four dressed correlators Gab can be found as 
the formal solution of the retarded Cauchy problem, with bare Greenians as initial data and the self-energies as the 
sources. Indeed, integrating the first term of each these equations twice by parts, we find for all four elements of Gab 



G^o{x,y) = [(^lfdY.^:^'^Gret{x,Ol''G^o{^.vh''Gadv{^l:V)~ I d^^dSGret{x,mA^-V)Gadv 

01 J 01 J 01 



('7,1/), 



(3.47) 



Goo(x,2/) = /rfS](«)dE('')G,et(x,e)7^Goo(^,77)7"Ga<i.('?,2/)- (d^idSGret{x,a±G^^ + ^^^{^.l^)]Gadv{^l,y). (3.48) 

11 J 11 J 00 

At first glance it is ambiguous to use the l ocal differential operator Gq^ in Eq.(3.48) without indication of direction 
it acts. Nevertheless, as follows from ( 3.42 ), 



Gret° oGadv 



Gr 



G 



(0) 



^Gadv — (1 + Gret O '^ret) O Gadv — Gret O (1 + '^adv ° Gadv) 



— —Gq + Gret O So O Gadv — 0, 

Dab and Mab 



and both directions lead to the same answer. 

The corresponding equations for boson fields are obtained by replacing Gab ^ Dab and M.ab ^ ^aBi and 

*d^dT,^ in Eqs. (jS^l through (|1|). 
The equations ( 3.47 ) and ( 3.48 ) as well as their copies for boson fields are the basic equations of relativistic quantum 
field kinetics. They are identical to the initial system of Schwinger-Dyson equations, but have the advantage that 
the time direction is explicitly singled out. There are two terms of different origin that contribute to any correlator 
(and, consequently, to any observable). The first term retains some memory about the initial data. The length of 
time for which this memory is kept depends on the retarded and advanced propagators. The second term describes 
the current dynamics of the system. A comparison of these two contribution allows one to judge if the system indeed 
has two scales. 

4. Perturbation theory for quantum field kinetics 



Any reasonable diagram technique should allow one to assemble certain subsequences of the bare diagrams into 
larger pieces, viz.^ the irreducible elements of the skeleton diagrams, such as exact Greenians, self-energies and vertices. 
The initial standard form of the Schwinger-Dyson equations (3.15), ( 3.16| ) and (3.22), ( 3.23| ) is an example of such 
selective summation. To derive these equations we have employed the functional approach, but any dressed field 
correlator, self-energy or vertex contains the infinite series of its perturbative expansion in powers of the coupling 
constant. We can reconstruct the formal perturbation series starting from the skeleton form of the exact Schwinger- 
Dyson equations. However, it is not expedient to proceed formally and we shall not attempt it here. For example, 
an explicit series of radiative corrections to the free asymptotic propagation is never retained in calculations. In the 
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course of renormalization this series is absorbed into the physical mass of the field. The theory of the deep inelastic 
e-p scattering that is based on the operator product expansion (OPE) selectively sums perturbation series in order to 
emphasize the dominance of light cone dynamics, and introduces its own irreducible elements, the structure functions. 
Just as the masses of free particles, they are taken from the data, and one should not try to obtain them directly 
from a perturbative expansion. Thus a physical {versus mathematical) study begins when we make the decision on 
how to sum the perturbation series or, equivalently, what skeleton elements are not to be expanded in a series. The 
final choice is the matter of taste. 



We wish now to perform the perturbative expansion of the skeleton equations (3.15), (3.18) and (3.22). Their matrix 
form accounts for a possible instability of the initial ground state with respect to interaction, but their mathematical 
form is precisely the same as in standard field theory where the initial ensemble corresponds to a pure and stable 
ground state. Any other structure would be surprising, as only the ordering of the operators, regardless of its type, 
is important. The presence of additional indices does not change the topology of the diagrams. 

Eqs. ( |3l5| ), ( pl8| ) and ( p^ , supplemented by the expressions for the self-energies, can be iterated in a standard 
way, giving rise to a formal perturbation series in the powers of coupling constant. Graphically, the series reproduces 
the Feynman diagram expansion with only one difference: each vertex acquires additional index. This kind of 
expansion proves useful when we solve local problems, such as photon or dilepton emission from a system with slowly 
varying macroscopic parameters. These parameters are built into the definition of the density matrix and should not 
change over the time of emission. 

We note, however, that this direct expansion is not most rewarding when the composite system undergoes a strong 
and short impact. Its parameters will change dramatically within a very short time, and th e me mory about the 
density matrix of the true initial state is obviously lost. It is then profitable to rearrange Eqs. ( ^l8| ) and ( |3^ ) (for 
example, "rotating" the basis of correlators, as in Eq. ( 3.30| )) in order to emphasize other details of the dynamical 
process. Since we are interested in processes on the light cone, the leading singularity of the retarded propagators is 
the main element which should be preserved in course of renormalization. 



4-1. Perturbation theory for the local emission problem 



The quantity which we shall choose to calculate in this section is the observable rate of the electromagnetic emission 
(real 7 or virtual 7*). It is expressed via the tensor of the electromagnetic polarization V^Ql—k) of the QGP, given 



by the Eq. (3.1!;). Let us expand it in a perturbation series up to order as. The polarization tensor 'P'^" contains 
three irreducible elements: two quark correlators and the electromagnetic vertex dressed by the strong interaction. 
They are given via their functional definitions, and obey the inhomogeneous integral equations ( p.l5| )-(3.21). The 
free terms of these equations correspond to the bare quark correlators, averaged over the ensemble represented by 
the density matrix (of non- interacting particles as a first approximation). The density matrix introduces occupation 
numbers which define the weights of initial and final states for elementary processes. 
The Born term emerges when all the irreducible elements are considered as bare: 



d'^p 



i-k) = ze'N^ J ^YGw{p~k)YGM- 



(4.1) 



To obtain corrections of order as to the electromagnetic polarization, one should begin by iterating Eqs. (3.15)-(3.21) 
to the same order. This iteration assumes that we write down equations (3.15) for the correlators with the bare 
Greenians and with self-energies computed to the first non- vanishing order. The vertex is considered to be bare. As a 
result of these approximations, we obtain the first group of radiative corrections. A second group appears if we keep 
field correlators bare, but include the first order corrections to the vertex. Eventually, such an expansion results in 
the following expression: 



V^i^i-k) = le^N, J A^7'^Gio(p - k)YGoi{p) + e'g 
d^pd^q 



c^F E (-1) 

R,S=0 



R+S 



-DsR{q)YGiR{p - k)^^Gm{p + q - kh'Gosip + q)lxGsi{p) 



where the correlators Gio and Gqi should be taken in the form 

Goi,10 — Goi,10 — GretMoi^ioGadv + GretMretGoi^iQ + Goi^ioMadyGadv 



(4.2) 



(4.3) 



All diagrams contributing to Vio{—k) up to order ae are depicted in Fig.|l|. The sum of these diagrams is the squared 
modulus of the coherent sum of amplitudes of the real processes of photon emission with first virtual corrections. These 
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amplitudes are given at Fig.g. The first process is the direct annihilation of a bare qq pair to a (virtual) photon, 
the next three diagrams are radiative contributions to the same process. Next, there follow four diagrams of qq- 
annihilation accompanied by gluon emission and absorption, and two diagrams for Compton scattering of a quark (or 
antiquark) and a gluon with photon emission. The first loop and the last four loops in Fig.^ are due to real processes, 
while the others are due to vertex and mass corrections. The correlators corresponding to initial or final states are 
marked by dashed lines. Lines with arrows correspond to retarded propagators, and arrows inside the loops show 
retarded self-energies. 

A complete calculation is a subject of separate paper. Here we shall only outline the main ideas and results. The 
technique of QFK allows one to perform all calculations without the assumption of thermal equilibrium. Not only 
can the statistical weights of the initial and final states of the real processes be of arbitrary form (though compatible 
with slowly varying macroscopic parameters), but we can also find all the virtual corrections consistently with the 
same quark-gluon background distributions. The general rule is that infrared divergences caused by the integration 
over the low momenta domain (of the Bloch-Nordsieck type) cancel out between real and virtual processes, despite 
the nontrivial population of the phase space. It is not even necessary to introduce any artificial intermediate cut-offs. 
Formally, for any given as-order the coUinear divergences survive, exactly because the phase space of the initial and 
final states is populated, and the balance required by the KLN theorem does not hold. However, the physical origin 
of these divergences allows one to conclude that they are shielded by the quark and gluon distributions: infinitely 
long free propagation in the dense system is impossible. 



4-3. Perturbation expansion for light cone processes 



As was already mentioned in Sections 1 and 2, QFK is equally applicable to the description of any inclusive 
processes, including those dominated by the light cone dynamics. These are deep inelastic e-p scattering, or deep 
inelastic collisions of two hadrons or nuclei. In this case the inclusive cross-section of DIS is expressed via the same 
electromagnetic polarization tensor, but the kinematic region is different. 

The cross-section for inclusive quark and gluon production, which might initiate evolution towards a QGP, is also 
expressed via the off-diagonal self-energies of the colliding system. For example, for gluon production, we have 



-^=J2Spp,„S^cHp,\,a)c{p,X,a)S = Y, j d^xd^^-^e^^h'^^^-rWlY^ {x,y)] (4.4) 



where the gluon polarization tensor Hqi is obtained by averaging with the density matrix of the system of the two 
colliding nuclei. The latter is the direct product of two independent density matrices for each nucleus. Thus, to first 
approximation, 

Klip) = ^9l{ j ^^S{k + q - p){Kr/(fc + q, -q, -k)T>l,Zd' (k)K'U-k - q, q, k)D[^,^f){q)] + (A ^ B)}. (4.5) 

where each of the gluon field correlators is averaged with the density matrix of one nucleus, A or B. The pre- 
collision dynamics of the field correlators T>l^^^f,{k) and D^^y,'^(g) is naturally described by equations similar to 



(3.46). The conditions for nuclear propagation along the light cone allows us to suggest that the retarded propagators 
contribute to the process mainly via their light cone singularity, and can be used without radiative corrections as a 
first approximation. In this case we may use the form 

-^01,10 — -^01,10 + -^01,10 '-'ret ^^01 so'-' adv ' V^-"^ 

where Dq'P^q describes the population of the gluon states used as initial data (at, for example, some scale Qq), and 

^oPw is the density of gluon states in the initially unpopulated continuum. 

The full description of the problem and the calculations is the subject of separate paper. Only the main ideas and 
results will be outlined here. For the sake of simplicity, let us consider the case of pure glue-dynamics. In that case, 



Eq. (4.5) yields the following evolution equation: 



n(;r(p) = -^9'r J ^ - ^) [Dretik)n,oik)Dad4kr'' V,^^^J,{-p,p - k, -k)D*';J,^{k - p) (4.7) 

If this equation is rewritten in light cone variables, and projected onto the specific measurement of the unpolarized 
e-p DIS experiment, it reproduces the well known Altarelli-Parisi evolution equation for the gluon structure function. 
If wi{p) is defined as the scalar coefficient in the decomposition 
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n>^'^ip)^gt^-^W,{p) + ^W2{p) + ... , (4.8) 



then the gluon structure function of the LLA reads 



where G{x,Qq) is the phenomenologic "initial" dist ribution at the scale Qq 



Here, the novel point is that equations like (4.7) are of the evolution type, even in their initial coordinate form. 
Further, they may be derived without reference to any particular process. The equations are of a ladder form, and 
the ordering of the ladder cells by the Feynman x and virtuality Q' is a consequence of the initial retarded ordering: 
lower X and higher Q' correspond to later times in the evolution. A parallel study of the two processes, p-e DIS, and 
p-p or A- A collisions, shows that one may use the new evolution equations to obtain quantities (like wj^(p)) which 
are common to these two processes. Further, one avoids intermediate parton language. 

7. Conclusions 

This paper presents a formalism which allows us to solve various problems of the QGP dynamics using the same 
technique. On the one hand, it allows one to calculate the rates of such processes as photon and dilepton emission, or 
heavy quarks production in the QGP. On the other hand, it is a good tool for the study of extremely non-equilibrium 
process, such as the first hard collision of two heavy ions. 

From the most general point of view, this technique is capable of replacing the standard Feynman approach, in any 
situation where the latter is applicable. In QFK, we begin with the density matrix generated by the Fock operators 
from the pure vacuum state of non-interacting fields. The reader can easily rewrite some chapters from a standard 
textbook on QED in terms of this approach as an exercise. In some points the new approach is physically more 
intuitive than the Feynman technique: at the tree level, retarded propagators appear naturally in every place where 
we usually put them in by hands (following common sense). The cutting rules and unitarity come to be a trivial 
consequence of the matrix structure: both S and S'^ contribute to any observable from the very beginning. 

For the case of true thermal equilibrium, this kinetic approach also reproduces all that is obtainable from the 
standard Matsubara formalism. All global relations like the FDT come to be a trivial consequence of thermal 
equilibrium and the matrix structure of the Schwinger-Dyson equations, but none of global relations are is used in 
the definition of the observables. For the Gibbs ensemble, QFK describes dynamical fluctuations. 

The essential new point is that the theory is not confined to a specific form of the one-particle distribution. The 
plasma may be far out of thermal and/or chemical equilibrium. Nevertheless we can consider both real and virtual 
processes in a self-consistent manner. Usually this consistency is lost in computer simulations of the emission from 
the non-equilibrium plasma j^]: only tree level matrix elements are taken into account. In the next paper we shall 
give an example: for dilepton emission from a non-equilibrium "plasma," a self-consistent account of the radiative 
corrections changes the answers from what is obtained using a naive cut-off of the infrared singularities. We have also 
found a significant difference in the case of an equilibrated plasma. 

The QFK approach allows one to formulate new principles for computing the distributions of quarks and giuons 
created in the first hard interaction of two heavy ions at high energies. It essentially employs an initial resummation 
of the perturbation series for the probabilities, and allows one to describe two different high energy processes, viz., 
e-p scattering and nuclear interactions, in the same terms. These processes are considered as two versions of the same 
phenomenon - deeply inelastic scattering of composite systems. The calculations can be performed without reference 
to parton phenomenology. 
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FIG. 1. Diagrams for the rate of dilepton production. The dashed line crosses the Greenians representing densities of the 
initial or final states. Arrows label the retarded and advanced propagators and show the latest time. 



FIG. 2. Processes participating in the dilepton production in the a^-order. Notation the same as in Fig.l. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9412362vl 




Fig. 2 



